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Abstract 

Direct  methanol  fuel  cells  (DMFCs)  have  significant  potential  to  become  a  leading  technology  for  energy  conversion  in  a  variety  of 
applications.  However,  problems,  such  as  methanol  crossover  reduce  the  efficiency  and  open  circuit  voltage  of  the  cells.  The  novel  design 
of  flowing  electrolyte-direct  methanol  fuel  cells  (FE-DMFCs)  addresses  this  issue.  Methanol  molecules  are  effectively  removed  from  the 
membrane  electrode  assembly  (MEA)  by  the  flowing  electrolyte,  and  the  unused  fuel  can  be  utilized  externally. 

In  this  paper,  a  general  3D  numerical  computational  fluid  dynamics  (CFD)  model  is  established  to  simulate  methanol  crossover  by 
convection-diffusion  in  the  FE-DMFC.  Illustrations  of  methanol  concentration  distribution  and  methanol  molar  flux  densities  are  presented, 
and  the  performance  is  compared  to  conventional  DMFCs.  The  results  indicate  that  methanol  crossover  can  be  reduced  significantly.  A 
parameter  study  is  performed  where  the  influences  of  anode  fuel  feed  concentration,  electrolyte  channel  thickness  and  electrolyte  volumetric 
flow  rate  on  methanol  crossover  are  evaluated.  In  addition,  effects  of  various  electrolyte  channel  orientations  are  determined.  According  to 
the  simulations,  counter  flow  is  the  superior  choice  of  channel  orientations  to  minimize  crossover. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 

Keywords:  Direct  methanol  fuel  cell  (DMFC);  Methanol  crossover;  Numerical  modeling;  Simulation;  Flowing  electrolyte;  Proton  exchange  membrane  (PEM) 


1.  Introduction 

The  research  efforts  and  interests  in  fuel  cell  technology 
and  development  are  increasing  rapidly.  Governments,  com¬ 
panies  and  universities  worldwide  are  gradually  adapting 
to  this  new  technology.  Due  to  high  efficiency  and  low 
emissions,  fuel  cells  display  a  great  variety  of  potential  appli¬ 
cations.  The  absence  of  moving  parts  and  an  extremely  simple 
mechanism  make  fuel  cells  very  competitive  for  small-scale 
applications.  For  such  applications,  low-temperature  fuel 
cells  with  proton  exchange  membranes  (PEM)  are  being 
adopted.  Currently,  efficient  hydrogen  PEM  fuel  cells  domi¬ 
nate  the  market,  but  the  problems  with  hydrogen  storage  and 
distribution  are  severe.  Liquid  fuels  are  much  easier  to  handle 
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although  the  conversion  in  fuel  cells  is  more  complex.  One 
of  the  best  available  liquid  fuels  in  terms  of  system  overall 
efficiency  (also  called  well  to  wheel  efficiency)  is  methanol 

[1]  that  can  be  utilized  directly  in  a  direct  methanol  fuel 
cell  (DMFC).  However,  the  DMFC  cannot  yet  compete  with 
hydrogen  fuel  cells  for  transportation  purposes,  mostly  due  to 
its  comparatively  low  cell  efficiency  and  expensive  catalysts 

[2] .  On  the  other  hand,  for  systems  where  size  and  weight 
are  more  important,  it  exhibits  high  potential.  The  power 
densities  and  energy  densities  of  DMFCs  are  superior,  even 
when  compared  to  newly  developed  lithium  ion  batteries  [2]. 
The  current  obstacles  are  low  cell  efficiency  and  high  cost, 
factors  that  must  be  improved  to  enable  commercialization. 

The  efficiency  of  direct  methanol  fuel  cells  is  reduced 
by  two  main  bottlenecks,  slow  anode  reaction  kinetics  and 
methanol  crossover.  Methanol  crossover  can  be  particularly 
severe  in  DMFCs  with  polymer  electrolytes.  Investigations 
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Nomenclature 

A  area  (m2) 

c  concentration  of  methanol  (molm-3) 

cv  specific  heat  at  constant  volume  ( J  kg- 1  K- 1 ) 
C2  concentration  of  methanol  at  the  cathode  side 
of  the  PEM  (molm-3) 

D  solute  diffusivity  (m2  s-1) 

A E  activation  energy  (J  mol-1) 

Fo  Fourier  number 

i  cell  current  density  (A  m-2) 

k\  constant  of  linearity 

kc  constant 

kp  constant 

ki  constant 

K  constant  related  to  effective  hydraulic  perme¬ 
ability  (m2  s-1  atm-1) 

Lpem  thickness  of  membrane  (m) 

Lx  length  of  MEA  in  FE  flow  direction  (m) 

Ly  thickness  of  FE  channel  (m) 

Lz  width  of  MEA  transverse  to  FE  flow  direction 

(m) 

h  methanol  molar  flux  density  (mol  s-1  m-2) 

P  pressure  (atm) 

q  internal  heat  source  (W  m-3) 

R  universal  gas  constant  (J  mol-1  K-1) 

S  source  of  species  (mol  s-1) 

t  time  (s) 

T  temperature  (K) 

v  fluid  velocity  (ms-1) 

V  volumetric  flow  rate  (m3  s-1) 

v  direction  of  FE  flow 

y  direction  across  FE  channel 

z  direction  of  FE  channel  width 


Greek  letters 

a  thermal  diffusivity  (m2  s-1) 

k  thermal  conductivity  (W  m- 1  K- 1 ) 

X  coefficient  of  electro-osmosis  (number  of 

methanol  molecules  dragged  by  each  proton) 
fj  chemical  potential  (J  mol- 1 ) 

fiv  kinematic  viscosity  (kgm-1  s-1) 

p  density  (kgm-3) 

Subscripts 

a  anode 

ave  average  value 

c  cathode 

FE  flowing  electrolyte 

max  maximum  value 

PEM  proton  exchange  membrane 


of  the  most  common  solid  polymer  electrolyte  material, 
Nafion®,  have  identified  that  methanol  can  penetrate  sulfonic 
acid  membranes  [3].  It  has  also  been  shown  that  methanol 
crossover  leads  to  a  significant  loss  in  oxygen  reduction  per¬ 
formance  at  the  cathode  due  to  a  mixed  potential.  Apart  from 
that,  crossover  is  a  substantial  loss  of  fuel  that  produces  no 
current  whatsoever.  Therefore,  overall  cell  potential  as  well 
as  fuel  efficiency  is  limited  by  this  effect.  Control  of  an¬ 
ode  fuel  feed  concentration  And  flow  rate  as  well  as  uti¬ 
lization  of  effective  methanol  oxidizing  anode  catalysts  are 
important  strategies  to  minimize  methanol  crossover.  There 
is  much  ongoing  work  in  this  field,  and  several  more  ad¬ 
vanced  methods  of  reducing  crossover  have  been  presented. 
For  instance,  the  electrolyte  material  can  be  replaced  [4,5] 
or  be  modified  chemically  [6,7]  to  achieve  lower  methanol 
permeation.  Shim  et  al.  [8]  suggest  the  use  of  a  methanol  im¬ 
permeable  barrier  between  the  electrodes,  allowing  protons 
to  pass  but  preventing  methanol  molecules  from  reaching 
the  cathode  catalyst  layer.  A  method  has  already  been  in¬ 
troduced  based  on  a  membrane  electrode  assembly  concept 
invented  by  Kordesch  et  al.  [9,10]  where  the  membrane  is  re¬ 
placed  by  a  flowing  electrolyte.  This  novel  design  efficiently 
removes  methanol  from  the  electrolyte  department  by  simple 
convection.  The  methanol  in  the  electrolyte  can  be  separated 
e.g.,  by  membrane  methods  and  consumed  or  a  small  amount 
may  be  used  up  in  a  loaded  control  sensing  cell  which  then 
automatically  adjusts  the  methanol  feed  to  the  system  [10]. 
Many  numerical  as  well  as  mathematical  models  of  the  con¬ 
ventional  DMFC  can  be  found  in  the  literature  [11,12],  but  to 
make  the  FE-DMFC  competitive  much  more  work  is  needed. 
A  mathematical  model  of  the  FE-DMFC  has  been  developed 
by  Golriz  et  al.  [13],  which  can  be  extended  by  support  of 
numerical  simulations. 

The  objective  of  this  paper  is  to  model  the  electrolyte  chan¬ 
nel  of  a  flowing  electrolyte  DMFC  and  numerically  simulate 
the  effect  of  methanol  crossover  for  different  operating  con¬ 
ditions.  Parameters,  such  as  electrolyte  flow  rate,  electrolyte 
channel  thickness,  fuel  concentration  and  channel  orienta¬ 
tion  will  be  studied  in  detail.  Simulations  of  the  conventional 
PEM-DMFC  will  also  be  performed  and  compared  to  the 
results  of  the  FE-DMFC  study.  The  model  is  intended  to  im¬ 
prove  the  knowledge  of  flowing  electrolyte  fuel  cell  perfor¬ 
mance  and  to  serve  as  a  vital  part  of  a  future  optimization 
study. 


2.  Flowing  electrolyte-direct  methanol  fuel  cells 

A  general  direct  methanol  fuel  cell  basically  consists  of 
a  membrane  electrode  assembly  (MEA)  supported  by  gas 
diffusion  layers  and  bipolar  plates.  The  MEA  contains  two 
porous  electrodes  (anode  and  cathode)  immersed  in  a  proton¬ 
conducting  electrolyte.  Fuel  (aqueous  methanol)  is  supplied 
to  the  anode  compartment  and  oxidizes  at  the  catalyst  layer 
to  form  carbon  dioxide,  protons  and  electrons.  The  protons 
diffuse  through  the  electrolyte  to  the  cathode  compartment, 
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Gas  diffusion  layer  Bipolar  plate 


Fig.  1 .  2D  cut  plane  of  a  typical  DMFC  module  with  proton  exchange  mem¬ 
brane. 

where  they  react  with  oxygen  to  produce  water  as  the  end 
product.  The  cell  reactions  are: 

CH3OH  +  H20  ->  C02  +  6H+  +  6e“  anode 
|02  +  6H+  +  6e_  — >►  3H2O  cathode 

CH3OH  +  |02  — >  2H2O  +  CO2  overall  cell  reaction 

The  electrons  generated  at  the  anode  flow  through  an  exter¬ 
nal  load  before  arriving  at  the  cathode,  thus  providing  useful 
electrical  power.  A  schematic  of  a  general  DMFC  with  proton 
exchange  membrane  is  presented  in  Fig.  1.  The  function  of 
the  membrane  is  to  physically  separate  the  electrodes,  transfer 
protons  from  anode  to  cathode  and  act  as  an  insulator  to  elec¬ 
tric  current.  A  good  electrolyte  shall  have  high  proton  con¬ 
ductivity,  low  methanol  permeation  and  low  electronic  con¬ 
ductivity,  and  shall  also  be  as  thin  as  possible.  Nation®  is  the 
commonly  used  material,  and  very  thin  films  can  be  produced. 
However,  the  problem  of  methanol  permeation  in  these  mem¬ 
branes  leads  to  reduced  cell  voltage  and  low  efficiency.  The 
new  concept  of  flowing  electrolyte- direct  methanol  fuel  cells 
(FE-DMFCs)  addresses  this  issue.  Here,  in  one  configura¬ 
tion,  an  electrolyte  channel  is  immersed  between  two  thin 
PEM  layers  as  in  Fig.  2.  The  flowing  electrolyte  (dilute  sul- 
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Fig.  2.  2D  cut  plane  of  the  corresponding  FE-DMFC  with  flowing  electrolyte 
(simplified  geometry,  not  to  scale). 


furic  acid)  effectively  removes  methanol  from  the  MEA  by 
convection,  and  hinders  its  arrival  at  the  cathode.  The  MEA 
of  the  second  generation  FE-DMFC,  with  a  50  mm  x  50  mm 
surface  area,  still  contains  two  thin  layers  of  PEM  (50  jam 
thick)  attached  to  the  electrodes.  However,  the  main  func¬ 
tion  of  these  layers  is  to  physically  separate  the  different 
channels  of  the  cell.  The  PEM-layers  could  be  considered 
optional  and  may  be  removed  or  replaced  by  a  thinner  and 
cheaper  material.  The  electrolyte  channel  width  is  a  target  for 
optimization.  Other  advantages  of  a  circulating  electrolyte 
system  are  simple  thermal  and  water  management;  the  possi¬ 
bility  to  remove  side  reaction  products  and  impurities,  such 
as  aldehyde  residues,  carboxylic  acids  or  other  intermediates 
produced  during  the  oxidation  of  methanol;  and  the  option 
to  adjust  the  electrolyte  flow  rate  relative  to  current  density 
or  fuel  concentration  to  minimize  fuel  crossover.  Moreover, 
methanol  transport  from  the  fuel  channel  to  the  anode  cata¬ 
lyst  layer  is  enhanced  by  the  flowing  electrolyte,  which  has 
been  shown  to  improve  performance  at  high  current  densities 
[14].  The  flowing  electrolyte  configuration  does  however  add 
some  extra  complexity  to  the  system,  e.g.,  electrolyte  piping 
and  pump. 


3.  Method 

As  mentioned  in  the  introduction,  the  objective  of  this 
work  is  to  create  a  numerical  model  of  the  electrolyte 
compartment  of  the  FE-DMFC  and  to  simulate  methanol 
crossover  for  various  operating  conditions.  All  parts  of  the 
fuel  cell  MEA  between  the  catalyst  layer  of  the  anode  and 
the  catalyst  layer  of  the  cathode  are  included  in  the  domain. 
The  catalyst  layers  of  the  electrodes  are  assumed  to  be  in¬ 
finitely  thin.  The  structure  of  the  electrolyte  department  in  a 
FE-DMFC  is  displayed  in  Fig.  2.  Notice  that  it  consists  of 
two  proton  exchange  membranes  with  a  flowing  electrolyte 
channel  in  between.  A  three-dimensional  model  is  necessary 
due  to  the  geometry  of  the  fuel  cell  in  general  and  the  elec¬ 
trolyte  channel  fluid  flow  in  particular.  For  this  purpose,  a 
computational  fluid  dynamics  (CFD)  model  is  an  obvious 
choice. 

Methanol  crossover  in  direct  methanol  fuel  cells  occurs 
due  to  the  driving  forces  of  chemical  potential  gradients, 
pressure  gradients  and  electro-osmosis  [15].  For  solutes,  the 
chemical  potential  gradient  equals  the  concentration  gradi¬ 
ent.  If  the  solute  diffusivity  is  assumed  to  be  independent  of 
concentration,  the  steady  state  methanol  molar  flux  density 
across  the  membrane  of  a  PEM-DMFC  is  described  by 

D  c2K  X 

^MeOH  =  - - Ac  -  - - A P  H - —i 

TpEM  Tpem  nt 

=  kc  Ac  +  kp  A P  +  kii  (1) 

According  to  the  model  of  a  conventional  DMFC  pre¬ 
sented  by  Wang  and  Wang  [12],  the  flux  due  to  the  pressure 
gradient  in  this  equation  is  small  compared  to  the  other  fluxes 
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and  can  therefore  be  neglected.  At  low  current  densities,  the 
diffusion  flux  dominates  and  electro-osmosis  flux  is  small. 
However,  at  high  current  densities  both  fluxes  are  significant. 
In  a  FE-DMFC,  with  a  fluid  flow  present  and  an  Eulerian 
formulation,  diffusion  is  replaced  by  convection-diffusion. 
The  process  is  more  complex  than  pure  diffusion,  but  the 
dominant  effect  is  a  suction  mechanism  that  enhances 
mass  transfer  by  convection-diffusion  through  the  anode 
side  membrane.  The  electro-osmotic  contribution,  on  the 
other  hand,  does  not  change  when  the  flowing  electrolyte 
convection  effect  is  introduced.  In  this  context,  especially 
if  we  restrict  the  operating  range  to  the  low  current  density 
regime,  the  contribution  of  electro-osmosis  may  also  be 
neglected.  The  flux  in  the  membrane  is  thus  simplified  to 

D 

^MeOH  =  ~~ - Ac  (2) 

TpEM 

In  this  paper,  only  the  convection-diffusion  contribution 
is  being  analyzed.  The  diffusivity  then  becomes  a  crucial 
parameter  for  the  quality  of  the  results. 

For  the  FE-DMFC,  the  temperature  dependence  of  the 
diffusivity  of  methanol  in  Nafion®,  AvieOH/PEM,  as  well  as 
methanol  in  water,  Dm&ou/u2o^  must  be  determined  (except 
for  its  high  proton  conductivity,  dilute  sulfuric  acid  is  as¬ 
sumed  to  have  the  same  properties  as  water).  In  general,  dif¬ 
fusivity  exhibits  exponential  temperature  dependence,  such 
as 


D(T)  =  Dref  exp 


Scott  et  al.  [16]  suggest  DmcOH/pem  =  4.9  x  10-10  m2  s-1 
at  333  K  and  DMeOH/H2o  =  2.8  x  10_9m2s_1  at  363  K 
as  reference  values  for  the  relationship  in  Eq.  (3),  and 
A£7R  =  2436K  for  Nafion®  117.  However,  several  other 
functions  can  be  found.  Kulikovsky  [17]  recommends  the 
reference  value  of  methanol  diffusivity  in  water  in  void  pores 
to  be  Dyq f  =  1 .58  10-9  m2  s-1  at  Tref  =  298  K,  and  use  the  tem¬ 
perature  dependence 


^MeOH/H2o(^)  =  A-ef  exp(0.026236(T  -  rref))  (4) 

This  is  in  agreement  with  the  reference  value  for  methanol 
in  water  given  by  Physics  Handbook  [18],  DMeOH/H2o  = 
1.4  x  10_9m2s_1  at  291  K,  although  both  differs  signifi¬ 
cantly  from  the  value  used  by  Scott  et  al.  [16].  For  methanol 
in  Nafion®,  Kulikovsky  [17]  uses 

^MeOH/PEM(^)  =  4.012  x  10-13  exp  (0.0243 12T)  (5) 

This  function  is  also  different  than  the  function  given  by 
Scott  et  al.  [16].  For  this  work,  the  diffusivities  given  by 
Scott  et  al.  [16]  are  chosen,  mainly  because  the  operating 
temperature  of  the  FE-DMFC  is  in  the  same  range  as  their 
reference  values.  Note  that  the  uncertainty  of  these  values 
can  be  considered  high. 

The  steady- state  velocity  profile  in  the  FE  channel  is  cal¬ 
culated  based  on  conservation  of  mass  and  momentum  in 


three  dimensions.  Conservation  of  mass  for  incompressible 
fluids  can  be  described  by  the  3D  continuity  equation 

V  •  v  =  0  (6) 

No- slip  boundary  conditions  and  the  Reynolds  number 
flow  imply  that  vy  =  0  and  vz  =  0  in  the  entire  channel.  Thus, 
assuming  that  the  flux  of  methanol  and  water  across  the  an¬ 
ode  side  membrane  is  small  compared  to  the  electrolyte  flow, 
fluid  flow  takes  place  in  the  v-direction  only.  Momentum  con¬ 
servation  is  described  by  the  3D  Navier-Stokes  equation  for 
incompressible  fluids 

dv  _ 

p—  +p(v-  V)v  =  -VP  +  pg  +  [iv  V2v  (7) 

at 

If  gravity  is  ignored,  the  steady-state  equation  can  be  writ¬ 
ten  as 

p(v  -  V)v  = —V P  +  jjvV2 v  (8) 

For  the  particular  geometry  of  the  FE  channel  (main  flow 
in  v-direction  from  inlet  to  outlet,  channel  width  much  larger 
left  to  right  (z)  than  between  the  top  and  bottom  membranes 
(y))  the  flow  in  the  middle  of  the  channel  can  be  modeled  as  a 
two-dimensional  steady-state  flow  between  two  plates  caused 
by  a  pressure  gradient  in  the  v-direction  (Poiseuille  flow).  The 
solution  of  Eq.  (8)  with  no-slip  boundary  conditions  is  a  linear 
pressure  gradient,  77,  along  the  electrolyte  flow  direction,  and 

vx(y)  =  P—y  {Ly  -  y )  (9) 

This  fully  developed  parabolic  velocity  profile  has  zero 
velocities  at  the  top  and  bottom  membrane  surfaces  and  a 
maximum  in  the  center  of  the  channel.  It  is  assumed  to  be 
representative  for  the  bulk  of  the  flow,  except  near  the  left 
and  right  channel  edges  at  z  =  0  and  Lz.  The  thickness  of  these 
edge  effects  is  assumed  to  be  of  the  same  order  of  magnitude 
as  the  channel  thickness  between  the  membranes.  The  veloc¬ 
ity  profile  is  assumed  to  have  a  parabolic  shape  where  edge 
effects  are  present,  with  zero  velocity  at  the  edge  increasing 
gradually  to  the  free  stream  velocity  as  one  moves  away  from 
the  edge.  The  free  stream  velocity  is  determined  by  Eq.  (9) 
independently  of  the  z-coordinate.  The  volumetric  flow  rate 
and  mean  velocity  are  obtained  by  surface  integration.  The 
Reynolds  number  is  calculated  by 

Re  =  PWL  (1Q) 

\Jv 

The  characteristic  length  of  the  FE  channel  is  its  thickness, 
Ly.  Reynolds  numbers  are  in  the  range  0.5  <Re  <  10,  thus 
indicating  laminar  flow. 

To  solve  the  mass  transfer  problem  in  the  presence  of  a 
velocity  field,  convection  as  well  as  diffusion  must  be  consid¬ 
ered.  Conservation  of  species  gives  the  convection-diffusion 
equation  for  time  dependent  solute  diffusion  in  3D  as 

—  +  v-Vii  +  V(-DVii)  =  S  (11) 

dt 
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For  a  dilute  species  in  a  liquid  solvent,  one  can  assume  that 
the  chemical  potential  equals  solute  concentration.  Thus,  with 
no  internal  chemical  reactions  or  sources  one  finally  obtains 
Fickian  advection-diffusion 

dc 

- 1-  v  •  Vc  +  V(— D  Vc)  =  0  (12) 

dt 


The  time  dependent  solute  convection-diffusion  equation 
in  the  presence  of  a  velocity  field,  Eq.  (12),  applied  to  the 
fuel  cell  MEA  domain,  is  solved  using  an  in-house  CFD 
software  called  SimManTec  Framework.  In  this  software 
code,  the  solute  convection-diffusion  equation  is  solved  us¬ 
ing  a  3D  transient  non-linear  space-time  method.  The  method 
treats  time  dependent  problems  like  stationary  problems  us¬ 
ing  space-time  meshes  that  discretize  space  and  time  in¬ 
dependently.  Usually,  finite  element  methods  use  a  semi¬ 
discrete  algorithm,  i.e.,  discretizes  space  using  FEM  and  use 
what  can  be  considered  a  finite  difference  method,  such  as 
Crank-Nicholson  in  time.  In  problems  where  the  domain  un¬ 
dergoes  deformations  in  time,  the  classical  finite  difference 
approach  may  become  difficult  to  use.  The  space-time  do¬ 
main  can  be  considered  a  four  dimensional  continuum,  where 
time  is  treated  as  the  fourth  dimension  space  variable.  But 
since  information  is  transferred  from  earlier  to  later  times, 
the  standard  stepwise  approach  from  finite  differences  is  pre¬ 
served.  A  continuous  Galerkin  method  in  space  and  a  discon¬ 
tinuous  method  in  time  are  used  for  the  discretization. 

The  methanol  crossover  model  of  the  FE-DMFC  is  incor¬ 
porated  as  a  time  dependent  problem  with  constant  velocity. 
Time  step  size  is  chosen  according  to  the  accuracy  criterion 
0.1  <  Fo  <  10,  based  on  the  dimensionless  Fourier  number 


Fo  = 


D  At 

A  y2 


(13) 


For  the  solute  convection-diffusion  problem,  Ay  repre¬ 
sents  the  mesh  element  size  in  the  anode-to-cathode  direc¬ 
tion  across  the  FE  channel.  The  operating  temperature  is  set 
to  60  °C.  Dirichlet  boundary  conditions  are  applied  to  the  an¬ 
ode  and  cathode  surfaces  facing  the  membranes.  Either  con¬ 
stant  or  linear  methanol  concentration  in  the  range  of  1-4  M 
is  applied  at  the  anode  surface,  and  zero  concentration  at  the 
cathode  surface.  Note  that  the  anode  surface  concentration  is 
expected  to  be  somewhat  lower  than  the  fuel  feed  concentra¬ 
tion,  due  to  the  concentration  drop  across  the  fuel  diffusion 
layer.  Also,  zero  concentration  at  the  cathode  indicates  that 
all  methanol  molecules  reaching  the  cathode  surface  are  be¬ 
ing  oxidized  instantly  at  the  catalyst  sites.  Neumann  bound¬ 
ary  conditions  with  zero  molar  flux  density  are  applied  to  all 
other  surfaces. 


4.  Results  and  discussion 

4.1.  Mesh  convergence  analysis  and  model  validation 

To  get  accurate  results  using  the  convection-diffusion 
solver  it  appears  to  be  crucial  to  employ  elements  that  are  not 


too  distorted  or  elongated.  If  the  aspect  ratio  of  the  element 
blocks  becomes  too  large,  the  solution  may  not  converge.  For 
the  geometry  of  the  MEA,  which  has  thin  membrane  films  of 
aspect  ratios  up  to  1000,  this  is  particularly  important.  When 
meshing  the  MEA,  this  implies  that  the  number  of  elements 
becomes  very  large.  For  accuracy  purposes  on  one  hand  and 
computational  time  on  the  other,  a  good  compromise  is  found 
with  eight  elements  of  eight-node  bricks  across  the  electrolyte 
channel.  The  error  using  this  mesh  is  expected  to  be  less  than 
0.2%  compared  to  an  ideal  mesh. 

The  model  is  validated  for  the  case  with  zero  FE  velocity, 
which  is  a  straightforward  analytical  calculation  using  Eq. 
(12).  With  zero  velocity,  the  steady- state  solution  is 


h  =  —D  Vc  =  —D 


Ac 

Ly 


(14) 


This  constant  methanol  molar  flux  implies  a  linear  concen¬ 
tration  distribution  across  the  membranes  and  the  stationary 
electrolyte.  The  gradient  is  higher  in  the  membrane  material 
(low  diffusivity)  than  in  the  flowing  electrolyte  (high  diffu- 
sivity)  and 


Cn  ~  Cc  =  Acpem  +  AcpE 


(15) 


In  this  equation,  Acpem  represents  the  sum  of  the  concen¬ 
tration  drops  across  both  membranes  in  the  MEA  of  the  FE- 
DMFC.  A  simulation  of  the  model  is  run  for  a  single  slice  in 
the  z-direction.  The  development  of  the  expected  concentra¬ 
tion  profile  with  time  is  monitored  and  steady-state  is  verified. 
Total  methanol  molar  flux  is  determined  by  post-processing 
surface  integration  applied  to  the  membrane  parts.  The  av¬ 
erage  molar  flux  across  the  anode  side  PEM  membrane  is 
computed  by 

^ave  —  ~  f  f  h  d A  (16) 

ApEM  J  J A 

This  gives  fiave  =  0.00671  mols-1  m~2,  which  is  constant 
throughout  all  parts  when  the  electrolyte  is  immobile.  The  de¬ 
viation  from  the  analytical  result  is  0.03%.  Also,  the  resulting 
concentration  distribution  is  clearly  linear  as  expected. 


4.2.  3D  simulations 


As  a  result  of  the  mesh  convergence  analysis,  3D  simula¬ 
tions  of  the  entire  domain  from  anode  to  cathode  with  cubic 
elements  become  computationally  expensive.  A  symmetry 
plane  in  the  middle  of  the  FE  channel  (z  =  25  mm)  is  identi¬ 
fied,  cutting  the  MEA  into  two  equivalent  parts  in  order  to 
reduce  the  number  of  elements.  Complete  3D  simulations  de¬ 
mand  a  trade-off  in  accuracy  but  for  illustrative  purposes  the 
precision  is  still  acceptable. 

Methanol  crossover  3D  simulations  of  the  MEA  part  of  the 
FE-DMFC  is  performed  for  the  case  with  electrolyte  flow 
rate  V  =  0.02  cm3  s-1,  FE  channel  thickness  Ly  =  0.6  mm 
and  constant  anode  concentration  ca  =  4  M.  The  time  step  size 
is  chosen  to  have  a  Fourier  number  in  the  range  0. 1  <  Fo  <  10, 
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and  a  time  dependent  simulation  is  run  until  steady  state  is 
reached.  Element  Courant  numbers  are  determined  by 

Co  =  — - -  (17) 

Ax 

For  the  time  step  size  chosen,  Co  =  20.  The  physical  ratio 
of  convection  to  diffusion  is  determined  by  the  Peclet  number 
for  the  problem 

Pe  =  LxVllve  =  3  x  104  (18) 

D 

Apparently,  convection  is  much  stronger  than  diffusion. 
The  results  obtained  consist  of  a  velocity  profile,  methanol 
concentration  distribution  and  molar  flux  density  at  all  nodes. 
As  described,  the  velocity  has  a  fully  developed  parabolic 
Poiseuille  profile  across  the  thickness  of  the  FE  channel.  In 
the  z-direction,  the  bulk  of  the  flow  is  uniform  except  near 
the  edges,  where  it  is  parabolic  due  to  edge  effects.  For  this 
velocity  profile,  the  volumetric  flow  rate  must  be  determined 
by  surface  integration  applied  to  the  channel  intersection 


The  total  volumetric  flow  rate  becomes  V  =  1.81  x 
10-8  m3  s-1,  which  is  slightly  lower  than  for  the  single  slice 
model  because  of  the  edge  effects  in  z.  Edge  effects  can  also 
be  seen  in  the  concentration  distribution.  Near  the  very  edge, 
where  the  velocity  approaches  zero,  the  concentration  is  lin¬ 
ear  across  the  FE  channel.  This  is  not  the  case  in  the  bulk  of  the 
fluid  flow,  where  the  flowing  electrolyte  efficiently  removes 
most  of  the  methanol  penetrating  the  anode  side  membrane. 
A  contour  plot  of  methanol  concentration  in  the  middle  of  the 
FE  channel  is  presented  in  Fig.  3.  The  figure  clearly  shows 
the  build-up  of  a  concentration  edge  effect  in  the  z-direction, 
which  appears  to  have  the  same  thickness  as  for  the  velocity. 
The  contours  indicate  that  the  bulk  of  the  flow  has  uniform 
concentration  distribution  along  the  z-axis.  In  the  case  with 
zero  fluid  velocity,  the  concentration  in  the  bulk  would  be 
equal  to  the  concentration  at  the  edge  (pure  diffusion).  Obvi- 
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Fig.  3.  Methanol  concentration  (shaded)  in  the  center  of  the  FE  channel 
(Ly  =  0.3  mm).  Fluid  flow  is  in  the  x-direction.  The  contours  represent  con¬ 
stant  concentration. 

ously,  convection  is  the  dominate  means  of  mass  transfer  in 
the  electrolyte  channel. 

Methanol  molar  flux  density  across  the  anode  and  cath¬ 
ode  surfaces  is  shown  in  Fig.  4.  Methanol  crossover  (cathode 
flux)  is  significantly  higher  near  the  outlet  of  the  channel  than 
close  to  the  inlet.  On  the  other  hand,  molar  flux  across  the 
anode  has  its  maximum  at  the  channel  inlet,  and  decreases 
along  the  flow  direction.  At  the  edge  (z  =  0),  the  flux  across 
the  anode  is  the  same  as  the  flux  across  the  cathode  (except 
near  the  inlet).  That  means,  where  the  fluid  velocity  is  zero 
no  species  are  removed  by  convection  and  the  only  transport 
mechanism  is  pure  diffusion.  Further  away  from  the  edge, 
convection  dominates.  At  a  short  distance  from  the  edge,  the 
flux  is  approximately  the  same  as  in  the  middle  of  the  channel 
(z  =  25  mm).  Note  that  the  anode  flux  is  significantly  larger 
than  for  pure  diffusion  due  to  the  suction  effect  at  the  anode 
surface.  The  suction  effect  may  enhance  mass  transfer  in  the 
catalyst  layer  of  the  anode  and  improve  reaction-rates  and  ef¬ 
ficiency.  Contrarily,  more  methanol  will  mix  into  the  flowing 
electrolyte.  It  is  thus  necessary  to  process  or  consume  this 
electrolyte  and  methanol  mixture  externally.  However,  it  is 
even  more  important  to  keep  methanol  away  from  the  cathode 
catalyst  layer.  For  this  purpose,  the  flowing  electrolyte  shows 
great  promise.  Except  very  near  the  z-edge  where  the  velocity 
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Fig.  4.  Anode  (a)  and  cathode  (b)  methanol  molar  flux  density  along  the  direction  of  fluid  flow  (x)  for  various  distances  from  the  edge  (z). 
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is  too  small,  methanol  crossover  is  reduced  to  a  minimum. 
The  average  molar  flux  density  across  anode  and  cathode 
is  calculated  by  surface  integration  as  in  Eq.  (17)  to  ha  = 
0.01284 mol  s-1  m-2  and  hc  =  0.001962 mol s-1  m-2.  In 
this  case,  85%  of  the  fuel  penetrating  the  anode  is  removed 
by  the  flowing  electrolyte,  which  is  more  than  three  times 
better  than  the  case  with  zero  electrolyte  velocity. 

As  a  comparison,  this  model  is  also  run  for  the  PEM- 
DMFC  case,  where  the  flowing  electrolyte  channel  is  re¬ 
moved  completely.  The  two  membranes  are  merged  to  one 
single  membrane  of  thickness  0.1  mm.  For  this  simple  ge¬ 
ometry,  methanol  flux  is  entirely  one-dimensional  and  all 
of  the  species  penetrating  the  anode  will  arrive  at  the  cath¬ 
ode.  The  results  for  ca  =  4  M  and  cc  =  0  M  are  na  =  hc  = 
0.0196 mol s-1  m-2.  The  methanol  crossover  in  the  FE- 
DMFC  is  reduced  by  90%  compared  to  this  result. 

4.3.  Parameter  study 

For  the  overall  performance  of  the  FE-DMFC,  a  3D  model 
of  at  least  half  the  MEA  is  essential.  However,  it  is  excessive 
for  a  detailed  parameter  study.  A  basically  two-dimensional 
model  of  one  single  slice  in  the  middle  of  the  channel 
(z  =  25  mm),  like  the  one  used  in  the  mesh  convergence  analy¬ 
sis,  is  instead  considered.  This  slice  with  only  one  element  in 
the  z-direction  is  assumed  to  be  representative  for  the  whole 
geometry.  From  a  system  point  of  view,  the  edge  effect  should 
be  included.  However,  since  the  bulk  of  the  geometry  dom¬ 
inates  the  performance  of  the  electrolyte  channel,  the  single 
slice  model  serves  the  purpose  of  a  parameter  study.  Simula¬ 
tions  on  the  single  slice  domain  are  carried  out  for  a  variety  of 
parameter  input  data.  Methanol  crossover  can  be  modeled  as 
a  function  of  several  variables,  such  as  operating  temperature, 
pressure,  FE  velocity  field  and  flow  rate,  FE  channel  dimen¬ 
sions,  methanol  concentration  at  anode  and  cathode,  current 
density,  cell  voltage,  etc.  In  this  preliminary  parameter  study, 
focus  is  on  three  of  these  variables;  methanol  concentration 
at  anode,  electrolyte  volumetric  flow  rate  and  FE  channel 
thickness.  The  temperature  is  constant  at  60  °C,  the  cathode 
concentration  is  assumed  to  be  zero  and  the  methanol  flux 
density  is  assumed  to  be  independent  of  the  other  variables. 
Thus, 

n  =  /(ca,  V,  Ly)  (20) 

Since  the  concentration  gradient  is  included  in  both  time- 
independent  terms  in  the  convection-diffusion  equation,  Eq. 
(13),  the  anode  concentration  is  expected  to  be  an  independent 
variable  and  is  therefore  treated  separately,  as  in 

n  =  /(Ca)  ■  f(V,  Ly)  (21) 

If  the  flow  rate  and  channel  thickness  are  held  constant, 
one  gets  a  linear  concentration  dependency 


(22) 
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Fig.  5 .  (a)  Average  methanol  molar  flux  density  across  anode  and  cathode  vs. 
anode  surface  concentration  for  FE-DMFC  and  PEM-DMFC.  (b)  Methanol 
molar  flux  density  fractions.  Conditions:  volumetric  flow  rate  0.02  cm3  s_1 
and  channel  thickness  0.6  mm. 


Plots  of  the  simulation  outcome  of  methanol  molar  flux 
density  across  anode  and  cathode  as  a  function  of  anode  con¬ 
centration  are  presented  in  Fig.  5a.  Methanol  crossover  in  the 
FE-DMFC  equals  cathode  flux.  Results  for  the  PEM-DMFC 
are  also  shown,  and  for  this  type,  methanol  flux  (crossover) 
is  the  same  at  both  electrodes.  Fig.  5b  displays  three  different 
flux  fractions;  crossover  fraction  for  FE-DMFC,  crossover 
fraction  compared  to  PEM-DMFC  and  anode  flux  fraction 
compared  to  PEM-DMFC.  These  plots  clearly  show  that 
methanol  molar  flux  is  directly  proportional  to  anode  surface 
concentration,  which  confirms  the  dependency  suggested  in 
Eq.  (22).  The  quantitative  results  are  similar  to  the  results  of 
the  complete  3D  simulation:  methanol  flux  across  the  cathode 
is  87.2%  lower  than  the  anode  flux,  and  methanol  crossover 
is  reduced  by  91.4%  compared  to  the  PEM-DMFC.  The  re¬ 
ductions  obtained  are  independent  of  anode  concentration. 
Additional  tests  are  performed  for  other  values  of  V  and  Ly, 
and  these  results  also  show  a  linear  dependency. 

The  results  obtained  for  ca  =  4  M  are  also  used  to  evaluate 
the  accuracy  of  the  single  slice  model.  Table  1  shows  average 
values  of  methanol  flux  from  the  single  slice  model  as  well 
as  the  full  3D  model,  which  is  used  as  the  reference  case. 


h  =  k\  •  ca 
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Table  1 


Comparison  of  the  results  obtained  with  the  single  slice  model  and  the  full 
3D  model 


Single  slice  model 

Full  3D  model 

Diff.  (%) 

ha  (mol  s-1  m-2) 

0.01314 

0.01284 

+2.3 

hc  (mol  s-1  m-2) 

0.001677 

0.001962 

-14.5 

V  (cm3  s-1) 

0.0197 

0.0181 

+8.8 

As  the  single  slice  model  ignores  the  edge  effects  at  the 
left  and  right  channel  edges,  it  is  expected  to  overestimate 
the  reduction  of  methanol  crossover.  This  is  confirmed  by  the 
14.5%  reduction  in  flux  across  the  cathode  compared  to  the 
3D  case.  In  the  same  manner,  anode  flux  is  slightly  higher. 
Moreover,  the  electrolyte  flow  rate  is  8.8%  higher  than  for 
the  full  3D  model  that  justifies  the  changes  in  methanol  flux. 
Considering  the  fact  that  the  models  operate  at  different  flow 
rates,  the  results  of  the  single  slice  model  are  quite  reasonable. 
In  practice,  edge  effects  can  be  avoided  by  constructing  the 
electrolyte  channel  a  few  millimeters  wider  than  the  catalyst 
surfaces  at  the  electrodes. 

Simulations  are  run  for  several  combinations  of  the  two 
other  variables,  FE  volumetric  flow  rate  and  FE  channel  thick¬ 
ness.  Four  different  levels  of  each  variable  are  used,  and 
methanol  concentration  at  the  anode  surface  is  constant.  El¬ 
ement  Courant  numbers  for  the  single  slice  model  are  in  the 
range  5  <  Co  <  50  and  the  problem  Peclet  numbers  are  be¬ 
tween  104  and  105.  The  results  of  these  tests  are  summa¬ 
rized  in  the  surface  plots  in  Fig.  6,  which  display  the  average 
methanol  molar  flux  density  as  a  function  of  both  variables. 
The  anode  molar  flux  increases  with  flow  rate  and  decreases 
with  channel  thickness.  The  methanol  suction  effect  appar¬ 
ently  is  enhanced  by  an  electrolyte  high  velocity  and  a  thin 
channel.  However,  the  appropriateness  of  inducing  methanol 
suction  is  a  matter  for  future  consideration:  it  may  be  prefer¬ 
able  to  avoid  methanol  penetration  of  the  anode,  thus  min¬ 
imizing  anode  flux  and  the  need  for  processing  of  unused 
fuel.  For  this  scenario,  a  wide  channel  with  a  low  flow  rate  is 
favorable.  Nevertheless,  anode  flux  in  the  FE-DMFC  is  still 
far  below  the  value  of  the  PEM-DMFC  (0.02  mol  s_1  m~2). 
Cathode  molar  flux  density  (methanol  crossover)  on  the  other 
hand,  is  most  certainly  a  means  of  minimization.  Fig.  7  dis¬ 
plays  the  achieved  reduction  of  crossover  compared  to  the 
equivalent  PEM-DMFC.  To  reduce  crossover,  high  FE  flow 
rates  and  wide  channels  are  to  be  used.  The  effect  of  in¬ 
creasing  either  one  of  them  appears  to  be  of  the  same  magni¬ 
tude.  Notice  that  except  for  the  cases  with  zero  flow  rates,  the 
crossover  reduction  is  more  than  65%  for  all  combinations 
being  tested. 

It  is  also  obvious  that  the  volumetric  flow  rate  and  chan¬ 
nel  thickness  variables  cannot  be  treated  independently.  The 
overall  output  effect  of  increasing  or  decreasing  one  of  the 
variables  may  be  the  same  independent  of  the  other  variable, 
but  the  magnitude  is  likely  to  vary.  For  an  overview  of  the 
effect  of  changing  channel  thickness,  Fig.  8  shows  results 
for  constant  flow  rate  V  =  0.02  cm3  s-1.  For  that  particular 
electrolyte  flow  rate,  methanol  crossover  is  kept  low  for  all 
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Fig.  6.  Surface  plots  of  average  methanol  molar  flux  density  across  (a)  the 
anode  surface  (line  AB  and  DC  are  atLy  =  1 .8  mm,  respectively);  (b)  the  cath¬ 
ode  surface  (crossover),  for  different  electrolyte  flow  rates  and  FE  channel 
thickness.  Anode  fuel  feed  concentration  is  kept  constant  at  4  M. 
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Fig.  7.  Reduction  of  methanol  crossover  compared  to  the  PEM-DMFC  for 
various  FE  channel  thicknesses  and  volumetric  flow  rates. 
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Fig.  8.  (a)  Average  methanol  molar  flux  density  vs.  FE  channel  thickness, 
(b)  Molar  flux  density  fractions.  FE  volumetric  flow  rate  0.02  cm3  s-1  and 
anode  surface  concentration  4  M  are  kept  constant. 

channel  thicknesses  analyzed.  For  Ly  =  0.3  mm,  the  crossover 
is  reduced  by  80%  and  for  Ly  =  0.6  mm  the  reduction  is  over 
90%  compared  to  the  PEM-DMFC.  The  plots  of  Fig.  9  also 
take  the  other  variable,  electrolyte  flow  rate,  into  account.  All 
four-channel  thicknesses  are  analyzed  and  compared  to  the 
PEM  case.  The  results  once  again  indicate  significant  removal 


Fig.  9.  Average  methanol  molar  flux  density  vs.  FE  volumetric  flow  rate  for 
different  channel  thickness.  Anode  surface  concentration  is  constant  at  4  M. 


Fig.  10.  Fractions  of  methanol  molar  flux  at  anode  and  cathode  vs.  electrolyte 
flow  rate  for  channel  thickness  0.3  and  0.6  mm.  Anode  surface  concentration 
is  4M. 

of  methanol  by  the  flowing  electrolyte.  Methanol  crossover 
is  reduced,  even  with  zero  flow  rates.  However,  there  is  a 
problem  with  the  insertion  of  the  electrolyte  channel.  There 
is  always  a  certain  resistance  in  the  electrolyte,  not  only  for 
methanol,  but  also  for  protons.  To  maintain  a  high  cell  volt¬ 
age,  the  rate  of  proton  transfer  from  anode  to  cathode  must  be 
kept  high  and  steady.  The  Ohmic  resistance  is  proportional 
to  channel  thickness,  and  it  is  important  to  avoid  unnecessary 
thick  channels.  Thick  channels  also  increase  system  size  and 
weight.  From  this  point  of  view,  it  is  most  likely  better  to 
increase  the  flow  rate  as  a  means  of  reducing  crossover,  and 
keep  the  channel  thickness  at  a  low  value. 

Fig.  10  displays  fractions  of  methanol  crossover  for  the 
0.3  and  0.6  mm  thick  channels  as  a  function  of  flow  rate.  For 
instance,  if  a  crossover  rate  less  than  20%  (80%  reduction) 
compared  to  the  PEM-DMFC  is  desired,  a  0.3  mm  wide  chan¬ 
nel  needs  a  flow  rate  of  0.02  cm3  s-1  whereas  for  a  0.6  mm 
channel,  0.01  cm3  s-1  is  sufficient.  For  crossover  rates  below 
10%,  these  flow  rates  should  be  doubled.  The  drawback  with 
an  increased  electrolyte  flow  rate  is  the  additional  pump  work 
needed.  The  pumping  power  must  be  monitored  and  included 
in  an  optimization  study  of  the  fuel  cell  performance. 

To  improve  the  performance  of  the  FE-DMFC,  the  flow¬ 
ing  electrolyte  channel  orientation  is  another  important  pa¬ 
rameter.  The  function  of  the  channel  itself  is  not  improved  by 
changing  its  orientation,  but  its  orientation  art  to  the  fuel  chan¬ 
nel  is  important.  The  anode  surface  concentration  obviously 
depends  on  fuel  channel  orientation.  Methanol  concentration 
in  the  fuel  feed  has  its  maximum  at  the  inlet  and  decreases 
gradually  towards  the  outlet.  Unlike  the  electrolyte  channel, 
the  fuel  channel  has  several  bends  and  a  complex  geometry. 
It  extends  over  the  anode  surface  like  a  maze.  For  the  model, 
the  fuel  channel  is  assumed  to  be  straight,  with  decreasing 
methanol  concentration  along  the  direction  of  the  flow.  This 
implies  a  linear  anode  surface  concentration.  There  are  three 
options  of  electrolyte  and  fuel  channel  orientations:  paral¬ 
lel  flow,  counter  flow  and  cross  flow.  Simulations  of  paral- 
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Fig.  1 1 .  Average  anode  (a)  and  cathode  (b)  methanol  molar  flux  for  the  cases 
with  constant  anode  concentration  3  M,  parallel  flow  and  counter  flow  with 
linear  anode  concentration  from  4  to  2  M.  Conditions:  FE  channel  thickness 
0.6  mm. 


Electrolyte  flow  rate  [cm3  s'1] 


Fig.  12.  Reduction  of  methanol  crossover  for  counter  flow  compared  to  par¬ 
allel  flow  and  constant  anode  concentration.  Conditions:  FE  channel  thick¬ 
ness  0.6  mm. 


lei  flow  and  counter  flow  are  performed  for  a  linear  anode 
concentration  starting  at  4M  and  ending  at  2M.  Simulation 
results  are  compared  to  the  case  of  constant  anode  concen¬ 
tration  ca  =  3M  in  Fig.  11.  Apparently,  methanol  crossover 
is  lower  for  counter  flow  than  for  parallel  flow  although  the 
anode  flux  is  enhanced.  For  parallel  flow,  the  crossover  rate 
is  even  higher  than  for  constant  anode  concentration.  Fig.  12 
presents  an  analysis  of  the  accomplished  crossover  reduction 
for  counter  flow  orientation.  The  results  indicate  excellent  im¬ 
provements  using  counter  flow.  For  an  electrolyte  flow  rate 
of  0.02  cm3  s_1 ,  the  reduction  of  methanol  crossover  is  26% 
and  15%  compared  to  parallel  flow  and  constant  anode  con¬ 
centration,  respectively. 


5.  Conclusions 

A  general  three-dimensional  numerical  CFD  model  of  the 
flowing  electrolyte-direct  methanol  fuel  cell  is  used  to  simu¬ 
late  methanol  crossover  by  convection  and  diffusion  for  a  va¬ 
riety  of  operating  conditions.  Methanol  concentration  distri¬ 
bution  and  methanol  molar  flux  densities  are  presented  for  the 
full  geometry  of  the  ME  A.  Compared  to  the  equivalent  pro¬ 
ton  exchange  membrane  direct  methanol  fuel  cell,  methanol 
crossover  is  greatly  reduced  due  to  enhanced  mass  transfer 
by  convection  (90%  reduction  for  reference  operating  con¬ 
ditions).  The  rate  of  crossover  is  shown  to  be  higher  in  the 
presence  of  concentration  edge  effects  near  the  outer  edges 
of  the  electrolyte  channel  than  in  the  bulk  of  the  electrolyte 
flow.  It  is  also  more  severe  near  the  outlet  than  at  the  inlet 
part  of  the  channel. 

For  a  thorough  study  of  various  parameters  and  operating 
conditions,  the  general  model  has  to  be  simplified  based  on 
a  slice  in  the  middle  of  the  electrolyte  channel,  which  ig¬ 
nores  the  edge  effects.  This  study  focuses  on  the  effect  of 
anode  surface  methanol  concentration,  flowing  electrolyte 
channel  thickness  and  volumetric  flow  rate.  The  average 
methanol  molar  flux  density  across  the  cathode  surface  (aver¬ 
age  methanol  crossover)  is  shown  to  be  directly  proportional 
to  the  anode  surface  concentration.  Combined  analyses  con¬ 
firm  that  methanol  crossover  is  inversely  proportional  to  both 
channel  thickness  and  flow  rate.  Thus,  wide  channels  and 
high  flow  rates  are  needed  to  minimize  it.  Due  to  the  am¬ 
plified  Ohmic  resistance  arising  for  wide  channels,  it  might 
be  a  better  strategy  to  increase  the  flow  rate  than  the  chan¬ 
nel  thickness.  But  with  increased  flow  rate,  it  is  important  to 
monitor  the  pump  work  needed  and  make  sure  it  does  not 
substantially  reduce  the  system  efficiency. 

The  convection  mechanism  also  affects  methanol  flux 
across  the  anode  surface.  High  flow  rates  and  low  channel 
thicknesses  lead  to  enhanced  anode  methanol  suction.  This 
improves  mass  transfer  in  the  anode  catalyst  layer,  which 
may  also  improve  the  reaction  rate.  However,  more  fuel  will 
be  mixed  in  the  electrolyte,  which  calls  for  external  process¬ 
ing  to  maintain  fuel  efficiency.  The  anode  flux  is  nevertheless 
significantly  lower  than  for  the  PEM-DMFC. 
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For  the  case  of  linearly  decreasing  anode  fuel  concentra¬ 
tion  from  inlet  to  outlet  of  the  fuel  channel,  different  elec¬ 
trolyte  channel  orientations  are  considered.  Parallel  flow  and 
counter  flow  are  both  compared  to  the  general  case  of  con¬ 
stant  anode  concentration.  To  reduce  fuel  crossover,  counter 
flow  is  shown  to  be  the  superior  channel  orientation. 

The  model  can  be  improved  further  by  adding  the  contri¬ 
butions  of  electro-osmosis  and  pressure  gradients  to  the  rate 
of  methanol  crossover,  which  is  an  ongoing  work.  More  ac¬ 
curate  values  of  the  diffusivity  of  methanol  in  Nafion®  and 
water  also  need  to  be  determined.  Further  work  includes  mod¬ 
eling  polarization  characteristics  and  determining  the  Ohmic 
resistance  of  the  electrolyte.  The  overall  objective  is  to  op¬ 
timize  the  flowing  electrolyte  channel  by  comparing  gain  in 
voltage  by  reduction  of  methanol  crossover  to  Ohmic  losses 
and  electrolyte  pump  work.  Additional  experiments  and  mea¬ 
surements  of  methanol  crossover  in  the  FE-DMFC  would  also 
be  desired. 
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